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ABSTRACT 

Knots or blobs observed in astrophysical jets are commonly interpreted as shock 
waves moving along them. Long time observations of the HST-1 knot inside the jet 
of the galaxy M87 have produced detailed multi-wavelength light curves. In this art¬ 
icle, we model t hese light curves using the semi-analytical approach developed by 
iMendoza et al.l (120091) . This model was developed to account for the light curves of 
working surfaces moving along relativistic jets. These working surfaces are generated 
by periodic oscillations of the injected flow velocity and mass ejection rates at the 
base of the jet. Using genetic algorithms to fit the parameters of the model, we are 
able to explain the outbursts observed in the light curves of the HST-1 knot with an 
accuracy greater than a 2-cr statistical confidence level. 
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1 INTRODUCTION 

T he jet in the g alaxy M87 was detected in the optical band 
by lCurtisl (119181 ). It is the closest Active Galaxy Nuclei with 
a redshift z = 0.004360 and has been extensively monitored 
in multi-frequency campaigns, particularly over the last dec¬ 
ade. Radio interferometry and high resolution optical and X- 
ray observations show the complex structur es formed inside 
the je t as close as ~ lOOpc from the nucleus ([Waters fc ZepJ 
120051 ). The most exotic of these structures, is a particular 
knot formed in 1999 and labelled HST-1. The evolution of 
HST-1 began to be closely followed in 2000 with the Chandra 
X-ray telescope (lHarris et al.ll2003l . 120061 . l2009h since it star¬ 
ted to develop a rapid increase on its X-ray emission, achiev¬ 
ing a maximum in 2005, corresponding to a factor of 50 as 
compared to the emission detected in 2000. After this max¬ 
imum, the emission decreases and is follo wed by t wo fu rther 
increments in 2006 and 2008. Ultraviolet (lMadridll2009l) and 
radio dChang et al.ll20icii ) observations show a similar beha¬ 
viour of its light curve. The whole emission of M87 presents 
an optical outburst in 2005 dMadridll2009l ) which is related 
to the maximum emission of the HST-1 knot in the same 
year. This strongly suggests that the outburst is produced 
by the strong emission of the knot. 

Knots in astrophysical jets are usually identified with 
internal shock waves travelling along the jet. These in¬ 
ternal shock waves can be produced by different mechan¬ 
isms: (a) intera ctions of the jet with an overdense medium , 
e.g. clouds (cf. iMendozal|2000I : IMendoza fe Lorigaitl[200ll) . 
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(b) bending of jets above a critical value (Mendoza l2000l : 
IMendoza fe Longairl 120021 1. and (c) Periodic variations of 
the injected velocity and mass at the bas e of the jet (e.g. 
Rees fe Meszarosl 1 19941 : I Jamil et al.l 120081 : IMendoza et aid 


20091 . and references therein). 


In the literature, the main contribution of the X-ray 
emission of the HST-1 knot is still under discussion and 
the interpretations v ary between an effect of a hot accretion 
disc with the corona dMarscher et al.ll2002l) and a particula r 
phenomena of a re-collimation shock ( Stawarz et al.ll2006l ). 
causing the impressive flare in X-rays. Later observations 
in radio revealed super luminal m otions in HST -1 being a 
well isolated knot from the nucleus dBiretta et al.ll 1 9991). dis¬ 
place d from the central engine by ^ 120 pc ( Cheung et al.l 
1200)1) . All this makes HST-1 the best studied knot for a 
possible internal shock mechanism inside a jet. It is also 
an ideal target to observe due to its proximity. The strong 
multi-wavelength emission from the jet and its knots allow 
to test the physics of knots and shock waves in the relativ¬ 
istic regime. 

Since relativistic outbursts are usually thought of as 
internal shock waves travelling along the jet, produced by 
periodic variations of the injected flow, it is quite natural 
to model the high emission light curve of the HST-1 knot 
as shock waves pr oduced by this mech anism. The semi- 
analytical model bv IMendoza et ahl (120091 ) (denoted as M09 
in what follows) has been quite useful in modelling not only 
outbursts associated to long gamma-ray bursts but also to 
the many outb ursts detected on the light curve of the blazar 
PKS 1510-089 dCabrera et al.ll2013l) . We show in this article, 
that such a model is also good for modelling and understand- 
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ing the multi-frequency features observed in the HST-1 knot 
of t he M87 galax y. 

lHarris et al.l (120091 ') found a quasi-periodic impulse sig¬ 
nature in the brightening and dimming of the core of M87. 
This was interpreted as a manifestation of past modulation 
of jet power, possibly by a local oscillation of the process 
that converts the bulk kinetic jet power to the internal en¬ 
ergy of the emitting plasma. This result reinforces the use 
of the M09 model in order to explain the formation and 
evolution of the HST-1 knot. 

The article is organised as follows. In section [2] we 
present the multi-wavelength observation campaigns of the 
HST-1 knot and its light curves features. In section [3] we 
present a b rief description of the hydrodynamical model de¬ 
veloped bv lMendoza et ahl (120091 ') and the system of dimen¬ 
sionless units in which it is useful to make comparisons with 
observations. The fi ts to the l i ght c urves using the hydro- 
dynamical model of lMendoza et all (l2009h are developed in 
the section [I] The result of our fits and a discussion of the 
obtained physical parameters of the model are presented in 
section [5] 


2 OBSERVATIONAL DATA 


The multi-frequency light curves were taken from three sep¬ 
arate datasets and are shown in Figure |T| X-ray observa¬ 
tions were taken from a multi-freque ncy program coord inat- 
ing Chandra and HST monitoring (lHarris et al.ll2009h . Ul¬ 
traviolet data are part of the s ame program and carried out 
during the years 1999 to 2006 dMadridll2009lh Finally the ra- 
dio data correspond s to observations with the VLBI at 2cm 
(Chang c t all 120 101. As mentioned in section [T| these all 
show a clear outburst with a maximum emission occurring 
in 2005, followed by a small outburst. After this, subsequent 
micro-outbursts differ from each other in the global decay 
of the light curve. Although all the observational data show 
the same morphology in the light curves, the spectral power 
law differs in each section of the spectrum from radio to X- 
rays, revealing that a si mple power law can not describe the 
whole spectra of HST- 1 (lHarris et al.ll2009h . 

We calculate t he flux in X - rays following the pro¬ 
cedure described by lHarris et al.l d2006l ~) and applied it to 
the observation al inte nsities of the HST-1 knot reported 
by lHarris et aid (120091 ). The flux in ultraviolet and radio 
wavelengths is calculated with a conversion factor between 
Jy to W/m 2 , using a reference wav elength of 225.55nm for 
the ultraviolet data (Madrid 2009j) and 2cm for the radio 
measurements dChang et aid 20101 '). 

We assume a negligible extinction factor and an iso¬ 
tropic emission of the source, located at a distance of 
1 6Mpc correspo nding to the distance to the galaxy M87 
d Jordan et al.l [20051 ') . With these assumptions we obtain a 
lower limit for the luminosity of the HST-1 knot in different 
wavelengths. 


3 MODEL 

Let us assume that periodic injections of velocity and mass 
flows are injected at the base of a ID relativistic flow mov¬ 
ing along a jet and consider a particular time on the ejec¬ 


S 

5 

>■4 

5 " 


1.02 
1 

0.98 
0.96 
0.94 
0.92 
0.9 
0.88 
0.86 
0.84 

1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 
Date(yrs) 



Figure 1. The figure shows multi-frequency luminosity curves of 
the HST-1 knot of the galaxy M87. All curves have been normal¬ 
ised to the maximum. From b ottom to top, the curves repr esent 
2cm radio (I Chang et ah Il2010h , 225 .5nm UV iMadridl (1200911 and 
Chandra 2Kev X-ray (lHarris et ah 1 120091) . 


tion process in which a fast parcel of flow is ejected after a 
slow one. A time later, the fast parcel will “overtake” the 
slow one and the flow will become multi-valued. In order to 
arrange this contradiction, nature creates an initial discon¬ 
tinuity of the hydrodynamical values which later develops 
into a working surface, i.e. a contact discontinuity bounded 
by two shock waves, that moves along the jet on the direc- 
tion of the flow as measu red from the central engine (see e.g. 
lLandau fe Lifshitzlll995h . 

The first ideas about radiative intern al shock waves in¬ 
side an astrophysical jet were develope d bv lRees fe Meszarosl 
ll994h : lDaigne fe Mochkovitchl (1 19981 1. Although several ex¬ 
tensions and particular aspects of the model have been 
presented in the l i terature (see e.g IP anaites cu et all Il999l : 
ISpada et al.l l200ll ; ISahavanathan fe Misra 2005ll . A semi- 
analytical description of this phenomenon was made by M09. 
These last model assumes periodic injections of mass and 
velocity at the base of the jet. Using mass and conserva¬ 
tions of the ejected material, it is possible to account for 
the kinetic power loss as the working surface travels along 
the jet, assuming that the radiation time scales are small 
compared to the characteristic dynamical times of the prob¬ 
lem. The pressure of the fluid is thus negligible and so 
the description of the flow can be well described by a bal¬ 
listic approximation. This assumption is valid if the flow 
within the jet is nearly adiabati c and non-turbulent (see 
e.g. ISahavanathan fe Misral I2005T ). In what follows we will 
use the model by M09 to describe the multi-wavelength light 
curve features of the HST-1 knot in M87. 

To follow the evolution of the working surfaces M09 
considered a source ejecting material in a preferred direction 
with a velocity v(t ) and a mass ejection rate m(r), both 
dependent on the time r as measured from the jet’s source. 
A further assumption is made such that the working surface 
is thin and mass losses within it are negligible. The energy 
loss E r by the working surface is given by E r = Eq — E ws , 
where Eo is the injected energy at the base of the jet and 
F ws is the energy inside the working surface. The kinetic 
power available within the working surface is then given by 
dE r /dt. If this power is converted efficiently into radiated 
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energy then the Luminosity L produced by the emission of 
the working surface is given by L = —dE r /dt. 

On the one hand, we assume that the injected velocity 
at the base of the jet is a periodic function of time, given 
by: 

v(r) = Vo + crj 2 sin (cur), (1) 

where the velocity Vo is the “background” average bulk ve¬ 
locity of the flow inside the jet and cu is the oscillation fre¬ 
quency of the injected velocity. The positive dimensionless 
parameter rj 2 measures the amplitude variations of the flow 
and is such that the oscillations of the flow are sufficiently 
small, in such a way that the total bulk velocity v(r) does 
not exceeds the velocity of light c. 

On the other hand, the mass ejection rate m injected 
at the base of the jet has the following periodic variation: 

m = rho + fi sin (fir). (2) 

where rho is the “background” average mass ejection rate 
and fi is the oscillation frequency of the mass ejection rate. 
The parameter fi is the amplitude of the injected oscillation. 

In the orig inal article by M09 and in further applic¬ 
ation s (see e.g. ICabrera et al.l 120131 : ICoronado fe Mendozal 
l2014l l the modelling of outbursts for long gamma-ray bursts, 
blazars and micro-quasars was performed under the assump¬ 
tion that fi = 0 and so m = const. Although this simplifies 
the number of free parameters of the model, it turns out 
that the light curve of the HST-1 knot in M87 cannot be 
modelled with such a simple assumption. 

In order to use the semi-analytical ballistic M09 
model on its more general form, we proceed as fol¬ 
lows. The model depends on six unknown parameters: 
Vo, r / 2 > rho, fi, and fi. To reduce the number of unknown 
parameters, we proceed as follows. 

The luminosity L depends on six dimensional paramet¬ 
ers: vo, cr/ 2 , w, rho, A, and fi. Additionally, the velocity 
of light c is an important dimensional parameter of the 
relativistic phenomena we are dealing with and so, it has 
to be added to the list of important dimensional quantit¬ 
ies of the problem. Since there are three fundamental inde¬ 
pendent dimensions, namely the dimensions of time, length 
and mass, Buckingham’s n-Theorem of dimensional analysis 
means that the luminosity can be described as follows: 

L — rn 0 c 2 L' (yo/c, rj 2 , (i/rn 0 , fi/w) . (3) 

In the previous equation, the dimensionless luminosity 
L' is a function of the four dimensionless quantities 
vo/c, rj 2 , js/rho, fi/w. In other words, the seven dimen¬ 
sional quantities for which the luminosity depends on, can 
be reduced to the problem of only four dimensionless quant¬ 
ities. 


4 FITS TO THE OBSERVATIONAL DATA 

The observed and theoretical luminosities, L abs and L t h re¬ 
spectively, can be fit to the observational data with the use 


of their dimensionless counterparts L' ohs and L[ h by rescal¬ 
ing them as follows. Both theoretical and observed dimen¬ 
sionless luminosities can be normalised to their maximum 
values. T 0 j 3S (T max ) and L, b ('T max ). i.e. 

Ubs := U (r' bS V Lth : = L' V (4) 

■^obs V ; obs,max / ■ L/ th\ / th ,max/ 

where the dimensionless times r/ bSjmax and 7th max corres¬ 
pond to the particular times where the observed or the¬ 
oretical luminosities reach a maximum value respectively. 
According to Buckingham’s n-Theorem of dimensional ana¬ 
lysis, the dimensionless time t' is related to the time r by 
the following relation: 

r = w _1 t . (5) 

In order to measure the observed and theoretical times in 
the same system of dimensionless units we normalised them 
to the time given by the FWHM of the outburst, i.e.: 

/ / 
rjp ._ _ ^"obs _ rjri _ _ ^~th _ /M 

° bs •“ rl(FWHM)’ th ^ r'(FWHM)' ( J 

The best fit of the theoretical luminosity L t h(T t h) 
to the observed light curve L 0 b s (T 0 b s ) yields a direct 
best value for the four dimensionless free parameters 
vo/c, rj 2 , ft/mo, fi/w. The quantity rho is obtained by us¬ 
ing © evaluated at one particular point of the light curve, 
which we choose as the point where the light curve reaches 
its maximum value. Once this last quantity is known, the 
value for the parameter (1 is hence inferred. The frequency 
uj is obtained using equation © evaluated at a particular 
time, which we choose as the time where the light curve 
reaches its maximum value. With this, the parameter fi is 
then inferred. 

The parameter calibration of the model is conceptual¬ 
ised as an optimisation problem and so, we propose to solve 
it using Genetic Algorithms (GAs), which are evolution¬ 
ary based stochastic search algorithms that mimics natural 
evolution. In this heuristic search technique, points in the 
search space are considered as individuals (solution candid¬ 
ates), which as a whole form a population. The particular 
fitness of an individual is a number, indicating its quality 
for the problem at hand. As in nature, GAs include a set 
of fundamental genetic operations that work on the geno¬ 
type, i.e. the solution candidate codificatio n, namely : muta - 
tion, recombination and selection operators iMitchehl (1 1998h . 
These algorithms operate with a population of individuals 
P(t) = x\, for a particular t iteration, where the fit¬ 

ness of each Xi individual is evaluated according to a set of 
objective functions fj(xi). This objectives functions allows 
to order from best to worst individuals of the population 
in a continuum of degrees of adaptation. Individuals with 
higher fitness, recombine their genotypes to form the gene 
pool of the next generation, in which random mutations are 
also introduced to produce a new variability. 

A fundamental advantage of GAs versus traditional 
methods is that GAs solve discrete, non-convex, discon¬ 
tinuous, and non-smooth problems successfully and so, 
they have been widely used in Ecolog y, Natural Resources 
Management, among other fields dLopez-Corona et al.l 


© 0000 RAS, MNRAS 000, 000-000 


















4 


Y. Coronado, 0. Lopez-Corona & S. Mendoza 



Date (yrs) 


Figure 2. The figure show the fits (lines) to the X-ra y data points 
of the light curve of the HST-1 knot observed by I Harris et al.l 
teoogh using the model by M09. The data points were divided into 
5 time sections marked by the dotted vertical lines corresponding 
to individual outbursts, and labelled XI, X2, X3, X4 and X5. The 
resulting calibration of the free parameters of the model by M09 
to the observed light curve are shown in Table [Q 

2Q13h with some ast rophysical applications (see e.g. 
Feigelson fe Babul 120121 ') . Our GA evaluated the luminos¬ 
ity function L t h (vo/c, rf, p/rho, fi/w) of M09 in order to 
compare numerical results from the model with the observed 
light curve L 0 b s using standard Residual Sum of Squares 
(RSS) as objective functions. All parameters were searched 
in the broadest possible range: 0.1 < vq/c < 0.999, 0.0001 < 
rf < 0.899, 0.001 < fi/rh 0 < 1.0 and 0.001 < fi/w < 20. The 
choice is consistent with the physical restriction of keeping 
subliminal the full bulk velocity of the flow v and to the fact 
that a large value of fi/iho would yield a huge unphysical 
luminosity value. A very large value of Q/uj produces large 
mass ejection oscillations, something not clearly visible from 
the light curves. This search parameter technique generates 
populations of 100 possible solutions over a maximum 5000 
generation search process, with a total of 500000 individu¬ 
als. The GA algori thms selected were: tournament selection 
with replacement (iGoldberg et alJll98Sl : l^stry _fc_Guklberg 
200lj), simulated binary crossov er (SBX) (|Peli_&^i<umai 
1993) and polynomial mutation (iDeb fe Kumadll995l : iDeb 
20011 1. The obtained final parameters were estimated by av¬ 
eraging the 500 best individuals. 

Direct inspection of the light curves in Figure [T] show 
that multiple outbursts occur during the per iod of observa¬ 
tion. A s such and following the procedure of ICabrera et alj 
(l2013ll we divided the light curves into individual outbursts. 
Two clear outbursts appear on all wavelength observations 
and an additional 3 mini-outburst were defined for the X-ray 
data -two before the main outburst and one at the end of the 
observations. The results of the GA explained above for each 
outburst are presented in Table [l] and the best fits to the 
light curves with these parameters are shown in Figures [50 


5 DISCUSSION 

Every modelling process goes an initial exploratory face in 
which a basic hypothesis space is set up. In this context, 
IWilliamsl (120141 ') found that a good modelling process should: 



Date (yrs) 

Figure 3. The figure show fits (lines) to the U V data points o f 
the light curve of the HST-1 knot observed by iMadridl (120091 ) 
using the model by M09. The data points were divided into 2 
time sections marked by the dotted vertical lines corresponding to 
individual outbursts, and labelled UV1 and UV2. The resulting 
calibration of the free parameters of the model by M09 to the 
observed light curve are shown in Tabled] 



Date (yrs) 


Figure 4. The figure show fits (lines) to the radio data points o f 
the light curve of the HST-1 knot observed by lChang ct all (1201(1 ) 
using the model by M09. The data points were divided into 2 
time sections marked by the dotted vertical lines corresponding 
to individual outbursts, and labelled R1 and R2. The resulting 
calibration of the free parameters of the model by M09 to the 
observed light curve are shown in Tabic [T] 

(a) stay as close as data as possible, (b) includes as much 
phenomenological information as possible and (c) keep as 
simple as possible. 

The parameter estimation of the model is quite close to 
the observational data (since it has a > 2<r confidence level 
value), with a simple ballistic model describing a complic¬ 
ated hydrodynamical phenomenon. 

At first sight, the curves seem not to properly adjust 
to many data points, as one should expect with such small 
observational uncertainties in the data. However, the time 
series represented by the light curve has many temporal 
gaps. Between these temporal gaps, the value of the inferred 
physical parameters may not stay the same, making the light 
curve to present mini-outbursts combined with different os¬ 
cillations. For example, the data points about 2007 in X-rays 
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ID 

Vo /c 

v 2 

7TIQ 

(lO^Moyr- 1 ) 

A 

(lO- 3 M 0 yr- 1 ) 

171 max 

(10- 3 M o yr- 1 ) 

ui Q 1 
(days) 

Q - 1 

(days) 

B min 

r max 

Tbulk 

Xi 

0.9631 

0.0360 

4.252 

0.415 

4.667 

53.8 

0.716 

2.67 

23.12 

3.72 

x 2 

0.9573 

0.0420 

3.476 

1.119 

4.596 

26.4 

4.361 

2.48 

27.32 

3.46 

x 3 

0.8156 

0.1839 

9.727 

3.943 

13.671 

374 

2322.5 

1.29 

31.32 

1.73 

X 4 

0.9713 

0.0275 

36.56 

3.351 

39.916 

118 

45.874 

3.03 

20.16 

4.2 

X 5 

0.9150 

0.0823 

4.511 

2.057 

6.569 

83.5 

36.0 

1.81 

13.49 

2.48 




(10- 6 M Q yr- 1 ) 

(lO- 6 M 0 yr— !) 

(10- 6 M o yr— !) 






UVi 

0.9020 

0.0976 

1.569 

0.588 

2.157 

296 

50.818 

1.68 

32.97 

2.32 

uv 2 

0.9762 

0.0228 

8.094 

3.338 

11.433 

16.1 

4.139 

3.31 

22.83 

4.61 




(10- 9 M Q yr- 1 ) 

(lO- 9 M 0 yr- 1 ) 

(10- 9 M o yr— !) 






Ri 

0.9450 

0.0536 

2.650 

2.033 

4.683 

195.5 

33.072 

2.21 

19.18 

3.06 

r 2 

0.9724 

0.0270 

9.156 

2.821 

11.978 

75.43 

31.21 

3.07 

27.12 

4.28 


Table 1 . Best parameter estimations using the X-ray, UV and radio light curves of the HST-1 knot of the galaxy M87. The resulting 
light curves are shown in Figures EH The fits were performed by dividing the light curves in time ID sections represented by the first 
column of the table. The quantity rh max corresponds to the maximum mass ejection rate discharged by the jet for a particular outburst. 
The values r imll . T max and IY u ik are the minimum, maximum and background (i.e. vq bulk “average” velocity of the flow) Lorentz 
factors of the flow. All parameters were obtained to a precision above a 2-a statistical confidence level. 


can be modelled as a series of mini-outbursts. But model¬ 
ling such a number of mini-bursts in a context of insufficient 
physical data represents an increase of unjustified additional 
hypothesis, despite the fact of an increment in sta tistical ac¬ 
curacy. As pointed out by iRoos fe Rakosl d2000l) it should 
be expected some sort of conflict between parsimony and 
realism. Nevertheless as models tends to incorporate more 
hypothesis, and become increasingly complex there is lose 
in transparency interpretation. 

Although equation m is dimensionally correct it 
doesn’t take into account the fact that an efficiency pro¬ 
portionality factor £ should appear in the right hand side of 
the equation, i.e. L = £moc 2 L' . This factor does not only 
depend on the ratio of the radiated luminosity to kinetic loss 
power inside the working surface, but also on the frequency 
of the emitted radiation. This is the reason as to why the 
inferred mass ejection rates for the same outburst differ so 
much at different wave-fengths. In other words, at best one 
should consider the values of the mass ejection rates in Table 
1 as lower limits. The inferred Lorentz factors for the bulk 
flow are ~ 1 — 4 reaching maximum values of up to ~ 30. 

The model by M09 has shown to be quite useful repro¬ 
ducing light curves of long gamma-ray bursts, blazars and 
micro-quasars. As we have shown in this article, the same 
model is also good for dealing with the light curve of the 
HST-1 knot of M87. Our modelling can be adjust more pre¬ 
cisely to the observed data by suitably performing more sub¬ 
divisions of the data set, essentially modelling many mini¬ 
outburst. Since no data is available for these mini-outbursts, 
their introduction would be speculative. In this sense, the 
current modelling can be i nterpreted as a baseline modelling 
dSchwab fe Starbuckl2013l ) that captures the key patterns in 
the empirical data and the associated physical processes. 
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